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Abstract 

We calculate explicitly the variation 5T C of the Bose-Einstein condensation 
temperature T c induced by weak repulsive two-body interactions to leading 
order in the interaction strength. As shown earlier by general arguments, 
5T C /T C is linear in the dimensionless product an 1 / 3 to leading order, where 
n is the density and a the scattering length. This result is non-perturbative, 
and a direct perturbative calculation of the amplitude is impossible due to 
infrared divergences familiar from the study of the superfluid helium lambda 
transition. Therefore we introduce here another standard expansion scheme, 
generalizing the initial model which depends on one complex field to one 
depending on ./V real fields, and calculating the temperature shift at leading 
order for large N. The result is explicit and finite. The reliability of the result 
depends on the relevance of the large N expansion to the situation N = 2, 
which can in principle be checked by systematic higher order calculations. 
The large N result agrees remarkably well with recent numerical simulations. 
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I. INTRODUCTION 



The effect of a weak repulsive two-body interaction on the transition temperature of a 
dilute gas Bose gas at fixed density has been controversial for a long time It has 

recently been established theoretically |7[] that T c increases linearly with the strength of the 
interaction parametrized in terms of the scattering length a. However the coefficient cannot 
be obtained from perturbation theory. In Ref. J?| a simple self-consistent approximation was 
used to derive an explicit estimate. 

In this article we first give a more direct derivation of the linear behavior, using general 
renormalization group arguments. Recognizing that the Hamiltonian of the system under 
study, which also describes the helium superfluid transition, is a particular example of the 
general N vector model, for iV =2, we generalize the problem to arbitrary N. This general- 
ization makes new tools available; in particular, the coefficient of AT C /T C can be calculated 
by carrying out an expansion in 1/N. Here we calculate explicitly to leading order in 1/N. 
The result happens to be independent of N, for non-trivial reasons. The calculation involves 
subtle technical points which are most easily dealt with by dimensional regularization. Sur- 
prisingly our result is in remarkable agreement with the most recent numerical simulations 

This paper is organized as follows: In Sec. II we lay out the basics of the problem. 
Then in Sec. Ill we present the general N vector model and analyze the behavior of the 
temperature shift by renormalization group arguments. Finally in Sec. IV we calculate the 
leading order contribution. 



II. CLASSICAL FIELD THEORY OF THE BOSE-EINSTEIN CONDENSATION 

We consider a system of identical bosons of mass m, at temperature T close to the critical 
temperature T c . The effective Hamiltonian of the system may be written as: 

U = f d 3 x (—Vi) ] (x) ■ Vif>(x) - fj,ifiUx)il)(x) + 2n?l a i) ] \x)i) ] '(x)tp(x)ip(x) \ , (1) 
J \2m m J 

where ip'(x) and ip(x) are the creation and the annihilation operators of the bosons, a is the 
scattering length, and \x is the chemical potential. Since we are interested in long wavelength 
phenomena, we have replaced the two-body potential by a delta function pseudopotential, 
with strength proportional to the scattering length a. We furthermore assume that a < A, 
where 



2nk 2 
\ rnksT 



(2) 



is the thermal wavelength. (In the following we use units h — ks — 1, and write simply 
A 2 = 2-n/mT.) 

To compute the effects of the interactions on the transition temperature, we write the 
particle number density as a sum of the single particle Green's function over Matsubara 
frequencies uj v = 2nii>T: 

n = -Tj^-^G(k^). (3) 
Above the transition, the single particle Green's function obeys the equation: 
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G 1 (k,z) = z + u 



2m 



E(k,z). 



(4) 



The Bose-Einstein condensation temperature is determined by the point where G 1 (0, 0) 
0, i.e., where £(0,0) = u. At this point, 



G~\k, z) = z-¥-- (E(/c, z) - E(0, 0)) . 
2m 



(5) 



At T c the spatial Fourier transform of the two-point correlation function at zero frequency 
diverges at zero momentum, and so does the correlation length. 
In the absence of interactions, u = at the transition, and 



n = C(3/2)/A* 



(6) 



where X 2 = 27i/mT®, and T c ° is the condensation temperature of the ideal gas. 

In the presence of weak interactions, the temperature of the Bose-Einstein condensation 
becomes the critical temperature of the interacting model, and is shifted by the interactions. 
From the theory of critical phenomena we know that the variation of the critical temperature 
in systems with dimension d below four depends primarily on contributions from the small 
momenta or large distance (which we refer to as the infrared, or IR) region. This property, 
which we later verify explicitly for d = 3, simplifies the problem, since to leading order the 
IR properties are only sensitive to the u v = component. 

As shown in Ref. J7|, in the dilute limit where only the u u = Matsubara frequency 
contributes, the shift in the transition temperature at fixed density, AT C = T c — T c °, can be 
related to the change An in the density at fixed T c by: 



2 An 

3 n 



where 



and 
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k 2 + U k 2 



U(k) =2m(E{k) - E(0)) . 



(7) 



(8) 



(9) 



Once restricted to their zero Matsubara frequency components, the fields if) and ip^ 
can be considered as classical fields, and the entire calculation can be cast in terms of 
classical field theory. It is then convenient to rescale the field if) in order to introduce more 
conventional field theory normalizations, and to parametrize it in terms of two real fields 

0i, (f>2- ip — VmT((f)i + i(j)2)- The partition function then reads 



Z = 

where S(<p) = H/T is given by: : 

1 



[d0(x)]exp [-5(0)], 



5(0) = 
where r = —2mT/j,, and: 



2 [d^(x)} + -r<f) 2 {x) 



u 
4! 



2 (x) 



dV 



(10) 



(ii) 
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u = 96n 2 — . (12) 
A 

In Eq. (|TT]) we have kept the dimension d of the spatial integration arbitrary in order to use 
dimensional regularization later. The single particle Green's function G(p) is related to the 
inverse two-point function T^(p) of the classical field theory by 

2mT 

G W = < 13 > 



The model described by the Euclidean action (|TT|) reduces to the ordinary 0(2) symmetric 
</> 4 field theory, which indeed describes the universal properties of the superfluid helium 
transition. As it stands this field theory suffers from UV divergences. These are absent 
in the original theory, the higher frequency modes providing a large momentum cutoff ~ 

VmT ~ 1/A. This cutoff may be restored when needed, e.g., by replacing the propagator 
by the regularized propagator: 

2mT 1 

fc2 * e k 2 /2mT _ I • V ) 

In fact, as we show later, since the shift of the critical temperature is dominated by long 
distance properties it is independent of the precise cutoff procedure. 

A second effect of the non-zero frequency modes is to renormalize the effective coefficients 
of the Euclidean action. This problem can be explored by writing the functional integral 
representation of the complete quantum theory and integrating over the non-zero modes 
perturbatively. The corrections generated are of higher order in a and can thus be neglected. 

Because the interactions are weak, one may imagine calculating the change in the transi- 
tion temperature by perturbation theory. However the perturbative expansion for a critical 
theory does not exist for any fixed dimension d < 4; IR divergences prevent a complete 
calculation. If one introduces an infrared cutoff k c to regulate the momentum integrals, one 
finds that perturbation theory breaks down when k c ~ a/A 2 , all terms being then of the 
same order of magnitude. To discuss this problem in more detail we now generalize the 
model to N component fields with an O(N) symmetric Hamiltonian. 



III. THE N- VECTOR MODEL. RENORMALIZATION GROUP 

We consider the O(N) symmetric generalization of the model corresponding to the Eu- 
clidean action (Tfl|). The field <p(x) then has N real components, and, e.g., 

N 

1=1 

The action S(<ft) is still given by Eq. (|TTD, now with an O(N) symmetry. The advantage of 
this generalization is that it provides us with a tool, the large N expansion, which allows us 
to calculate at the critical point |T0|-|T2"||. 

The goal is to obtain the leading order non-trivial contribution at criticality (in the 
massless theory) to 

N 

n = 2mT ^ (cj)f] = 2mT N p, (16) 
i=i 
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with 

- f ddk 1 
P ~ J (2n) d rW{k)' { } 

Here Sij/T^(k) is the connected two-point correlation function. 

We now recover, by a simple renormalization group analysis, the result of |7] that the 
change in p due to the interaction is linear in the coupling constant. We introduce a large 

momentum cutoff A oc y/mT ~ 1/A, and the dimensionless coupling constant g 

S = A d " 4 uoc (£f~ 2 . (18) 

At T c the inverse two-point function in momentum space satisfies the renormalization group 
equation || 

( A Jv + m §~g ~ Vi9) ) ^ A ' 9) = ° • (19) 

This equation with dimensional analysis implies that the two-point function has the general 
form 

T^(p,A,g)=p 2 Z(g)F(p/A(g)); (20) 
on dimensional grounds A.(g) is proportional to A, with 

m^-=m, (21.) 
m ?mm = . l . ( 2ib) 

Since /3(g) = -(4- d)g + (N + 8)g 2 /48it 2 + 0(g 3 ), (3(g) is of order g for for small g in d < 4; 
similarly rj(g) = (N + 2)^ 2 /(72(8vr 2 ) 2 ) + 0(g 3 ). Therefore 

Z(g)=expJ Q 9 ^dg' = l + 0(g 2 ); (22) 

to leading order Z(g) =1. The function A(^) is then obtained by integrating Eq. ( |21b|) , 



Kg) = 9 1/(4 - d) ^ 



exp 



9 A > < 1 ' 



K (3(g') (A-d)g' / 



(23) 



The scale A(g) plays a specific role in the analysis as the crossover separating a universal 
long-distance regime, where 

T^(p) oc p 2 -' r > p«A(g), (24) 

governed by the non-trivial zero, g*, of the /3-function, from a universal short distance regime 
governed by the Gaussian fixed point, g = 0, where 

T {2) (p)ocp 2 A(#)<p<A. (25) 
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However such a regime exists only if A(g) -C A, i.e., if there is an intermediate scale between 
the IR and the microscopic scales; otherwise only the IR behavior can be observed. In a 
generic situation g is of order unity, and thus A(g) is of order A, and the universal large 
momentum region is absent. Instead A(g) <C A implies 



P(g') (4-dV 



<1. (26) 



Since g (equal to a/A, see Eq. (0)) is <C 1, this condition is satisfied in the present situation. 

We now show that with this condition, AT C oc A(g). First, from the g = limit we infer 
that F(oo) = 1 (see Eq. (0).) In three dimensions the function F(p) behaves for large p as 

F(p) = l + 0(\np/p 2 ), (27) 

as can be seen directly from 3-d perturbation theory. Therefore the first correction to the 
density flTTD is convergent at large momentum and independent of the cutoff procedure, 

r d 3 p 1 / 1 \ 

Sp= J(^^{Fwm- 1 )- (28) 

Similarly the IR scaling result (Eq. [24]) implies that this integral is IR convergent. Setting 
p = A(g)k, we then find the general form 

6p = M9)[^lJ^-A; (29) 



(2n) 3 k 2 \F(k) 

the g dependence is entirely contained in A(g). For g small we conclude 

^ oc g oc an 1 ' 3 , (30) 

in agreement with J7[. It is important to note that both the perturbative large momentum 
region and the non-perturbative IR region contribute to the integral in Eq. (|29f) . Therefore 
we cannot evaluate it from a perturbative calculation of the function F(p). However, we 
now show, we can calculate 5p exactly in an 1/N expansion. 



IV. THE LARGE N EXPANSION AT ORDER 1/JV 

Critical phenomena can be studied in any fixed dimension by the now standard technique 
of large N expansion, where the large limit is taken at Nu fixed. To leading order the 
critical two-point function has simply its free field form. However a non-trivial correction is 
generated at order 1/N; one finds the inverse two-point function |13|.|9"|j. 

(2) 2 2 r d d q 1 ( \ 1\ / 1 \ 

[P)-P + ( 2 Ti) d (6/Nu) + B(q) \{p + qf ~ ) + \N* ) ' ( ' 

where B(q) is the one-loop contribution to the four-point function 

, [ d d k 1 

= J wm^f ■ < 32) 
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Note that in the large N limit, the chemical potential \x is proportional to 1/N. Using the 
large iV value of g*, and setting e = 4 — d, we can write B(q) for q — > 0: 

5(g) = &(e)<r - j^A- £ + 0(1/ A 2 ), (33) 



with 



6(e) = ~ • A/ 9 mi d/ n ^' ( 34 ) 
sin(7ra/2j 1 (a — 1) 

where A^ is the usual loop factor 

2 1 

^ = P^P)' iV3 = 2^- (35) 

For d = 3, 6(1) = 1/8. 

In the large iV limit the /3-function takes the simple form (3(g) = —eg(l—g/g*). Therefore 
the leading cutoff-dependent correction to B(q) combines with 6/Nu to yield (6/N)A~ £ (g), 
as expected from renormalization group arguments. This cutoff dependent correction can 
be neglected for small g. 

We evaluate 

_2_ r d d p 1 d d q ( 1 _ 1_\ 

P A^i (2tt) 2 V (6/#u) + 6( £ )g- £ ^(p + g) 2 g 2 J 1 j 

by keeping the dimension d generic and using dimensional regularization [[14]]. The integral 
over p is 



d d p 1 1 1 r(3 - d/2)T(d/2 - l)r(d/2 - 2) d 



(2vr) d p 4 (p + g) 2 (47r) d / 2 r(d - 3) 



Q 



-6 



1 r(d/2 - 1) 



7T 



g d ~ 6 . (37) 



(47r) d / 2 r(d-3) sin7rrf/2 



Note that the singularity at d = 3, which would apparently entail the vanishing of the 
integral, is cancelled in the subsequent q integral, which reduces to 

(2vr) rf (6/iVu) + 6(e)g- £ 4-d sin (vr(rf - 2)/(4 - d)) VAW ' V ; 

In the d = 3 limit the two integrations in Eq. (RH) yield (1/32tt 2 )(Nu/6). As expected, 
5p oc u: 

5p = -Ku, K = -^, (39) 

967T 2 

or in terms of the original parameters, 

Sp = ~. (40) 

It is instructive to repeat the calculation directly at d — 3. To do this, we first eliminate 
the constant term (which does not contribute), and write U(k), which equals r^ 2 ^(A;) — A; 2 , 
as 
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U(k) = -2N 



d 3 g 



B(q) 



(2vr) 3 1 + (Nu/6)B(q) \(k - q)< 



(41) 



From Eq. 



we find 



U(k) 



-327r 2 Nk— 



dx 



A 4 Jo kx + 2n 2 aN/X 2 \2 



*ln 



1 + x 



1 — X 



(42) 



Calculating the change in the density in leading order in U ~ 1/iV, one gets: 



dfc 



dx 



7T 2 A 2 JO JO 1 + Txfc \2 



x 



In 



1 + x 



1 — X 



- 1 



(43) 



where r = A 2 /(27r 2 aiV). The trick now is to exchange the orders of the k and x integrals. 
In three dimensions, however, the integrals are not absolutely convergent, preventing this 
interchange. Thus we introduce a regularization, with by inserting a factor k e in the k 
integral, and take the limit e — > + . With this factor we may exchange the orders of 
integration. For small e the k integral becomes 



oo he-l 
dk- 



1 + rxk e(rx) e 
The factor r e goes to unity, and the remaining x integral becomes 



(44) 



(45) 



For e = this integral vanishes identically. Thus we may replace x e by x e — 1 which goes 
to — elnx as e —>■ 0. The remaining integral is 





1 + X 




/d««-(fln 


-0 




1 — X 





dx In x ( — In 



1 + x 



X 



- 1 



7T 



(46) 



The factors of e cancel out, and we find, as before (Eq. (f4"U|)), 5p = —a/A 2 . 

Using this result in Eq. (|7|), we finally obtain the change in the transition temperature: 



8tt 



3C(3/2) A 
8tt 

3C(3/2) 4 /3 



an 



1/3 



2.33 an 1/3 , 



(47) 



Note that although the final result does not depend on N and therefore replacing iV by 
two is easy, the result is only valid for iV large. The result ( f4T|) is in remarkable agreement 

with the (N = 2) value AT C /T C ° w (2.2 ± 0.2)an 1 / 3 in the recent numerical simulations of 
Holzmann and Krauth [HI. 



V. CONCLUSION 



In this paper we have shown that the properties of the weakly interacting Bose gas 
remain dominated by the UV fixed point of the renormalization group equations up to very 
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large length scales; this is why we can still refer to the Bose-Einstein condensation when 
discussing the phase transition of the dilute interacting Bose gas. Renormalization group 
arguments also enabled us to confirm directly that the shift of the transition temperature at 
fixed density is proportional to the dimensionless combination an 1 / 3 for weak interactions. 
This result is non-perturbative, and the proportionality coefficient cannot be obtained from 
perturbation theory. We have therefore introduced a non-perturbative method, the large 
N expansion, which allows a systematic calculation of this coefficient as a power series in 
l/N, where eventually one has to set N = 2. Finally we have calculated explicitly the 
leading order contribution. The first correction is formally of order l/N multiplied by a 
function of the product aN which is fixed in the large N limit. Because the final result in 
three dimensions is linear in a, the l/N factor somewhat surprisingly cancels, and the result 
is independent of N. The value found is in remarkable agreement with the most recent 
numerical estimates. Whether this agreement is just a coincidence or reflects the smallness 
of the next order (of order l/N) correction can only be checked by an explicit calculation. 
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